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Abstract 


We develop a new model for flexible workforce management in environments with uncertainty in the demand for 
labor. In particular we model a workforce comprised of regular workers who have fixed schedules and may work 
overtime, and contingent workers whose working hours are flexible as specified by a contract over a finite planning 
horizon. Our model can represent a variety of contracts, including: temporary workers, on-call workers with guaranteed 
minimum pay, and comp-time arrangements. We formulate the model as an optimization problem that determines the 
regular and contingent worker pool sizes that minimize expected labor and backlog costs. Embedded within this 
problem is a dynamic programming problem of making optimal operational staffing decisions with respect to the 
utilization of contingent and overtime resources. Numerical examples demonstrate the effect that the timing of infor- 
mation has on the benefits of flexibility. We also derive structural results that may be exploited to reduce the com- 
putational effort required to use the model. 
© 2002 Elsevier Science B.V. All rights reserved. 
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1. Introduction and literature review 


Between 1990 and 2000 revenue in the temporary help industry grew from 17.0 to 63.6 billion dollars 
(not inflation adjusted), and annual employment increased from .99 million to 2.54 million (American 
Staffing Association, 2001). This rapid growth is but one sign of the dramatic changes the US labor market 
has been undergoing in the past decade. While the US labor market has traditionally been a leader in labor 
flexibility the same trends are appearing in Europe as well. For example, in the year 2000, 60% of the sales 
of Manpower, the leading staffing services provider, came from its European operations (Manpower Inc., 
2000). Throughout the economy, especially in services, firms rely more and more upon external contract 
workers to fill key positions that were once considered the exclusive purview of full-time permanent workers 
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Nomenclature 

Parameters 

V number of periods in planning horizon 

G guaranteed fraction of V periods contingent workers are paid for 

Crws Cows Cot COSt per shift of regular workers, contingent workers and overtime work, respectively 
Cb cost of a unit of backlog or unsatisfied demand per period 

CB cost of a unit of backlog or unsatisfied demand in the final period 

Cr fixed cost per contingent worker per planning horizon 


Tw, Tew, Hot Productivity of regular workers, contingent workers, and overtime work in units of work 
per shift, respectively 

OT,,,, OT, parameters determining upper limit of overtime shifts contingent and regular workers 
may work in a period, respectively 

OT(S) upper bound on the total overtime shifts in a period for staff S = (N,M) 


Decision variables 

N,M — number of regular and contingent workers staffed for the planning horizon 
O; the number of overtime shifts utilized in even stages ¢ 

Uy the number of contingent worker shifts utilized in odd stages t¢ 


Random variables 
n,, N actual and expected number of regular workers present in stage ¢ 


X; work in system at start of stage ¢, the net productivity of staff already committed, in units of 
regular worker productivity in a shift 

Ky unused portion of the total contingent worker guarantee at end of stage t— 1 

~; the amount of new work arriving to the system in stage ¢ for odd ¢, in units of regular worker 
productivity in a shift 

0, the amount of new work arriving to the system in stage ¢ for even ¢, in units of regular worker 


productivity in a shift 


(Applebaum, 1987; Clinton, 1997; Brotherton, 1995). At the same time a variety of flexible staffing ar- 
rangements have become more prevalent within firms. Some of these arrangements include part-time, 
comp-time ' and on-call work (Christensen, 1995). The causes of these changes are quite varied and 
complex and, being recent phenomena, are still under investigation. In fact the US Bureau of Labor 
Statistics only began collecting data on contingent workers in 1995 (US Bureau of Labor Statistics, 1995). 
When the adoption of these different internal and external staffing arrangements is driven by a firm’s 
uncertainty in its demand for labor we refer to them as contingent work arrangements, using the definition 
in Polivka (1989): “Any job in which an individual does not have an explicit or implicit contract for long- 
term employment or one in which the minimum hours worked can vary in a non-systematic manner’. 
As a result of the above changes in the labor market, managers now have many more choices in crafting 
a Staffing strategy that flexibly matches labor supply to demand. In this paper we formulate and analyze an 
operational model of a firm’s use of contingent labor that provides managers with a tool to determine how 
many full-time and contingent workers to hire. Furthermore we provide a framework for quantifying the 
effects of having more or less demand information when making short-term staffing adjustments. We be- 


' Comp-time is the practice by which workers are compensated for overtime in one week by shorter hours in the following week. 
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lieve such a tool is important for two major reasons. First, there is much anecdotal evidence that firms use 
contingent workers and find that they do not reduce costs (Nollen and Axel, 1995). This suggests either that 
they do not know how to use contingent workers effectively, or that they are not able to accurately evaluate 
the benefits of doing so. Second, from the authors’ experiences working with the United States postal 
service (USPS) and a financial services firm it is apparent to us that without the correct information systems 
in place to track labor demand, it is difficult to take advantage of the flexibility of contingent workers. Our 
numerical experiments reveal three interesting insights. First, we see that information and flexibility 
strengthen each other’s cost reduction effects. Second, we see that information and capacity act as com- 
plements, in the economic sense, in that increasing the demand information available increases the staffing 
levels selected. Thirdly, we see that increased contingent labor flexibility does not always decrease regular 
worker staffing levels. 

We formulate a model of contingent labor management based upon our observations of the back room 
operations of various services. We consider a workplace in which there is uncertainty in the demand for 
labor resources (on a daily or weekly basis), there is uncertainty in the supply of labor (absenteeism), work 
can be backlogged, and no inventory of processed work is allowed. We assume that the employer may draw 
upon three sources of labor: fixed, over-time, and contingent, to minimize costs of labor and backlogged 
work. This description could apply to many clerical operations (mail sorting, insurance claims processing 
etc.), delivery services, and repair services among others. The key characteristic is that work is created by a 
customer request. Explicit to our model is a representation of how demand information becomes available 
to the decision maker and a model of contingent labor contracts. 

A single-period model of the use of contingent workers appears in Abraham (1986) assuming an unlimited 
supply of flexible labor to a firm and no particular structure to the relationship. Berman and Larson (1994) 
model the use of contingent workers in an environment where demand is realized each day before decisions on 
the use of contingent workers must be made, there is no bound on overtime each day and no backlogging of 
unfinished work. A hierarchical model for planning the use of long-term, medium-term, and short-term 
flexibility is presented in Wild and Schneeweiss (1993). That paper models the use of short-timing (cutting back 
full-time hours), temporary workers, overtime, and crosstrained workers in an environment with uncertain 
demand for labor and absenteeism. It assumes that there is no backlogging or inventories and therefore daily 
staffing decisions are independent of one another. The model demonstrates the importance of the interplay 
between long-term planning decisions and the availability of resources for short-term staffing decisions. The 
model does not, however, represent different types of contingent worker contracts and it does not represent the 
way demand information timing affects their utilization. Our paper addresses both of these issues. 

There has been recent interest in modeling contracts between suppliers and buyers within the inventory 
control and supply chain management literature (see Donohue, 2000; Anupindi and Bassok, 1998, 1997; 
Chen et al., 2001 for example). While these models have shed some light on some of the important issues of 
supplier—buyer relations, they are not structured to be directly applied to the workforce management setting. 
An exception is Milner and Pinker (2001) which considers labor supply contracts but not at the tactical 
decisionmaking level done here. This particular combination of overtime, processing capacity constraints, 
guaranteed work, and lack of inventories has not been addressed in the broader production literature. 

In Sections 2 and 3 we, respectively, formulate and analyze our model. In Section 4 we present a set of 
illustrative numerical results. We conclude in Section 5. For the reader’s convenience a glossary of notation 
used in the model is also given. 


2. Model formulation 


We study a firm that must set staffing levels for a planning horizon of V periods. The firm decides how 
many regular workers N to hire and how many contingent workers © to contract from a contingent labor 
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supplier. Regular workers are each guaranteed a wage (including benefits) of c,j per period. The cost of the 
contingent labor is established by the contract. The contract specifies how many contingent workers M are 
contracted, a fixed payment of cp per worker, a wage rate (including benefits) of c., per period worked by 
each worker and a guaranteed number of person-periods of work MGV, where G is (0, 1]. 

In this section we formulate a mathematical model of the firm’s staffing problem that includes the use 
of regular workers, contingent workers, and overtime and allows for uncertainty in demand for labor, 
backlogging of unfinished work, absenteeism, differing productivity rates across classes of workers, and 
different demand information scenarios. To properly assess the actual labor and backlog costs incurred by a 
particular solution, we model the operation of the workplace with optimal labor resource allocation de- 
cisions made on the period level. We model this level of activity with a dynamic programming model in 
which the manager has two decisions each period: how many contingent workers to use and how much 
overtime to use. Over the entire planning horizon these decisions determine the total expected labor and 
backlog costs for the system. In each period new work arrives to the system and is combined with the 
existing backlog. The manager then draws upon the three labor sources to process the work. The regular 
workers are automatically used, while contingent worker and overtime usage are both decided by the 
manager. Unprocessed work in a period is backlogged to the next period. Since the manager may have to 
decide how many contingent workers to call upon before deciding how much overtime to use the manager 
may have more information about the period’s demand when making the latter decision. 


2.1. Definitions and assumptions 


Order of events. We define a time period to be an arbitrary, problem specific, unit of time. As discussed in 
Section 1, the model developed in this paper is geared toward periods of duration of the order of days or 
weeks, and within each period two distinct staffing decisions are made, how many contingent workers to use 
and how much overtime work to use. We assume the following sequence of events during each period: 


A first batch of work arrives and regular worker absenteeism is realized. 
The contingent worker usage decision is made. 

A second batch of work arrives. 

The overtime decision is made. 


By “work arrival” we mean demand information rather than physical arrival of work. Some initial 
information about the quantity of new work for a period may arrive before the contingent worker decision 
is made, while additional information arrives after the decision is made. For example at the end of one work 
day (Monday) a manager may know how many regular workers have asked for the following day (Tuesday) 
off and have a partial forecast of how much work there will be the next day (Tuesday). At this point she can 
decide how many contingent workers to call in advance to come to work the next day (Tuesday). During 
Tuesday the manager will see the true extent of the workload and can decide to use overtime. The shorter 
the notice the manager has to provide contingent workers the more demand information available to her 
when she makes that decision. To reflect the fact that each of the staffing decisions may be made with 
differing amounts of demand information and absenteeism information we split each period 7 into a first or 
odd stage ¢ = 2i— 1 and a second or even stage t+ 1 = 27, wherei=1,...,V with t=1,...,2V. 

For an odd stage ¢ one decides how many contingent workers to use, u, = u,(n,, @,), where n, (the number 
of regular workers present) and @, (the first arrival of work in the period) are realized quantities. For an 
even stage ¢ one decides how much overtime work to use, w, = «,(0,), where 0, (the second arrival of work 
in the period) is the realized quantity. 

Note that it is arbitrary to require the absenteeism of regular workers to be known before contingent 
worker decisions are made. The model could be formulated, with few changes, with the absenteeism in- 
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formation becoming available after the contingent worker decision is made (see Pinker, 1996). The alter- 
native formulation leads to similar results. 

Demand process. We assume that the ¢,’s are independent, that all the 0,’s are independent as well, and 
that the ¢,’s and 0,’s are independent of each other. To simplify notation in this paper we assume that the 
o,’s are identically distributed as are the 0,’s, although the model could accommodate non-homogenous 
demand processes. The sum of the two random variables ¢,_, and 0, defines the total exogenous demand for 
labor in the period spanning the stages ¢— 1 and 1, i.e., period t/2. 

Absenteeism. We assume that the distribution of the number of regular workers present in a period is the 
same in all periods. We also assume that regular workers are paid for periods they are absent. This is similar 
to practice, in that workers are given an allowance of sick and personal days. We also assume that there is 
no absenteeism among contingent workers. This means that the contingent labor supplier is compelled to 
always provide up to M contingent workers by the contract. 

Overtime. We consider any hours attributed to contingent workers after u, has been reached, 1.e., in even 
stages, as overtime and not counted toward the guarantee. For even stages t we decide w,, the number of 
overtime shifts utilized in stage t. Overtime each period is constrained such that 


0, < OT wW + OTM = OT(S), 


where WN is the expected number of regular workers who are present in a period, and OT,, and OT.y are 
parameters, in units of shifts, that determine the maximum allowable overtime an individual worker of each 
type may perform in one period. Since the bound on overtime is expressed in terms of the expected number 
of regular workers present and not the actual number, there can be cases in which individuals exceed the 
overtime bound. Therefore, the overtime bound is an approximation of practice, made in the interest of 
tractability. We assume that contingent and regular workers both have the same productivity, z.,, when 
performing overtime, also in the interest of tractability. 

Cost parameters. We use cp to denote the fixed cost for each contingent worker that is a member of the 
workforce for the planning period. It includes fixed component of compensation and fixed costs per worker 
for the firm, e.g., human resource departments costs, services available to all employees regardless of status, 
etc. Let c, be the per-shift cost of a regular worker that combines the benefits and salary. Let cey be the per- 
shift cost of a contingent worker that combines the prorated benefits and salary. Let c., be the cost per- 
worker per shift of overtime worked. Let c, be the penalty per period incurred by the firm for every unit of 
work backlogged in the first V — 1 periods, and let cg be the corresponding penalty for the final period of 
the planning horizon. We allow for a different backlog penalty in the final period to provide flexibility to the 
user of this model. 


2.2. Pool sizing problem formulation 


We now formulate the pool sizing problem as 
(P) min C(S) = CrwNV + cpM + CovMGV + Ey, n, {fi(0, MGV) } 
s.t. S=(N,M) integer and B(S) > (0,0), 


where B(S) is a set of side constraints that may restrict the proportion of workers that may be of one type or 
the other. The cost of hiring N regular workers is c,ywNV, and c;M + cyyMGV is the up front cost of M 
contingent workers. The expression Ey, »,{fi(0, MGV)} is the expected cost of making optimal tactical 
staffing decisions over the planning horizon given a staffing level S = (N,M), starting with 0 units of work 
in the system, and MGV guaranteed shifts of contingent workers unused in stage 1. The state of the system 
in stage ¢ is defined by x, units of work (net productivity of committed staff) and x, guaranteed shifts of 
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contingent labor. In the following we use the notation (y)* and [y]* interchangeably to denote the positive 
part of y. To simplify the notation we define a new random variable ¢) = $, — twn;, which is new work 
arrival in odd stages net work done by the regular workers who are present. ¢, can take on negative values. 
We can then define f;(x,, «,) recursively as follows. 

The cost-to-go function in stage 2V is: 


Sov (Xo, Kay) = min {Cot@2v + cp(X2v — Tot@2v) }, 
Ve 


Woy < min[x2y /7.,OT(S)], where OT(S) = OT wN + OTM. 
For even f: 


Sil, Kr) = min {CoO + Co (%1 — Horr) + Ey, esi (%es1, Ker) ]}, 
ow, < min|[x,/7,,OT(S)], where OT(S) = OTN + OTM. 
In even stages ¢ the cost-to-go is a combination of overtime charges c,,, backlog penalties c, and expected 
future costs of operating the system. 
For odd t:  f,(x;, K+) = min { Cow [Ur —K)* + Eos Feri (Xe415 Kis1)]}5 O0<u,<M. 
In odd stages ¢ the cost-to-go when u, contingent workers are solicited is determined by the number of 


remaining guaranteed shifts plus the expected future costs of operating the system. 
From even stage t to odd stage f+ 1 state transitions occur as follows: 


/ 
Xp = Xp — Mor Ms + Pint and Ki41 = Ky. 


From odd stage ¢ to even stage f+ 1 state transitions occur as follows: 
Xt = [xy = Towly + Geea ls 


Kil = [Ky = ul". 


We have defined a two-dimensional state space (see Fig. 1) (x;,«,), where x, is the unused portion of the 
total contingent worker guarantee, MGV, at the end of stage ¢ — 1, and x, is the work in the system at the 
start of stage ¢, net the productivity of staff already committed for stage t. Note x, may take on negative 
values for odd ¢. If x; is positive it means that the backlogged work in the system exceeds the regular staff 
processing capacity. If x, is negative it means that there are more regular workers present than needed for 
the backlog carried over from stage t — | and therefore the staff can be assigned to work on the exogenous 
work arriving to the system in the second stage. 


KA 
MGV Oo OR 
Remaining contingent 
guarantee t ae! 
T Even stages 
Odd stage Pte tal 
> 
X 
(0,0) Work in System : 


Fig. 1. State space for tactical problem dynamic program. 
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3. Analysis 


Solving the staffing problem requires the determination of a cost minimizing staffing level S = (NV, M), 
where the evaluation of a particular solution S involves the optimal solution of a stochastic dynamic 
programming problem. In this section we present results that simplify the computational effort necessary 
for the dynamic optimization. The key results are summarized as follows: 


The cost-to-go functions of the dynamic optimization are convex in the state variables. 
Finding the optimal overtime and contingent worker utilization policies are convex minimization prob- 
lems in each stage. 

e There are easily computed bounds on the optimal u, and w, in each period. 


The key to efficiently optimizing the model is to take advantage of the fact that the size of the dynamic 
program embedded within the model is directly proportional to the size of the contingent labor pool M and 
the guarantee G being evaluated. Therefore we prefer to evaluate the DP for smaller values of MG. The 
results we derive in this paper allow us to calculate bounds on the decision variables M, u,, and @, using 
only calculations in which the « dimension of the DP state space is eliminated. 


3.1. Dynamic optimization 


Proposition 1. f,(x,,«,) is convex, in (x,,,), for all t, and the choice of u, and w, are convex minimization 
problems for all t. 


Proof. See Appendix A. 


3.1.1. Overtime policy 
For even ¢ we can rewrite the cost-to-go function as 


Co 
Silt, Kr) = Lt, Kr) + Xn 
ot 
where 
L(x, ky) = min {T(by, K:)}, 


[x;-OT(S) 01] * <b <M 


ot 


T (by, Kr) > (0 ee “ + Ey, fivi(br TF Prats K-41) 


by = Xp — Torr. 


The change of variables means that selecting the optimal amount of overtime to use is equivalent to se- 
lecting an optimal amount of backlog, b,, to accept or tolerate in the system. Let us denote the backlog that 
minimizes I’;(b,,«,) for a particular x, in period ¢ to be $,(x,) which exists due to the convexity of f,(x;, K;). 
That is, £,(x,) is the minimizer of I,(b,, «,) without regard to the overtime usage constraint in period ¢ or the 
work in the system x,. Once we have determined f,(x,), the optimal overtime w*(x,,«,) for a particular state 
is found by comparing f,(«,) with x, and [x, — t..OT(S)]” as follows: 
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Lemma 1. Given [,(«,) and x, > 0 the optimal overtime usage is: w* = (x; — b*)/(tot), where 


Xt if B,( Kr) 2 Xt, 
Di (x1, K1) = 4 Bie) if = TOT (S)|" <BACK) <1, 
[x; — MrOT(S)]* if B,(«,) < [x, — tox OT(S)]*. 


Proof. Follows from Proposition 1 and discussion above. 


When cy > Cot/Mot, then B,(x,) = 0 for all «, because the per period cost of a unit of backlog is greater 
than the overtime cost of eliminating it. Using the above characterization of the optimal overtime policy we 
can prove that the following bound on @, holds. 


Proposition 2. When cy < Cot/Tor; 


@; S (%1 — B21, 0))/Tor- 


Proof. See Appendix A. 


3.1.2. Contingent worker policy 

We note that the size of the state space of the model increases with MGV. In particular when G = 0 we 
have a one-dimensional problem that is solvable relatively quickly. We use the optimal contingent worker 
utilization policy from this special case to compute lower bounds on contingent worker utilization in 
general. We first state the following preliminary result regarding the optimal w,. 


Proposition 3. When x, = 0 the optimal contingent worker usage in odd stages t, U,, is characterized as fol- 
lows: 


0 if x, <X;, 
U,(x;) = if X,+ Mio =x, > X, 
M if x >X,+Mrw, 


where X, is the minimizer of (—Cow/Tew)X: + Eo,,, fis ((Xr + 6141], 0). 


Proof. See Appendix A. 


The X, defined in Proposition 3 can be computed using standard backwards recursion. If «, > 0, then 
contingent workers are cheaper to use than if x, = 0, and therefore, u*(x) > U,(x) for all values of x, 1.e., 
the U,(x) are an easily computed lower bound on w(x) in each odd stage. 


3.2. Pool-sizing problem 


In our experience with minimizing the objective function of problem (P), C(S), we have found that it is 
convex and therefore a gradient search or binary search procedure always converges. However, we have not 
been able to prove this mathematically. Intuitively one would expect that the cost reduction of adding 
additional staff on the dynamic programming component of the objective function diminishes as S in- 
creases. 
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We also make the following observations without proof: 

(Ol) The optimal regular worker staff size N* should be bounded above by the optimal N when no 
contingent workers are available, N*(M = 0). 

(O2) The following condition: N + tyM + 2.OT(S) > E[¢, + 0;] should be true to prevent instability. 

(O3) It should be the case that M*(N) < M*(N — 1), where M*(N) is the optimal value of M given a fixed 
N. 

The implication of (O3) is that by searching through the space of feasible S starting with large 
values of N and small values of M we can improve the lower bound on M*(N) with solutions to smaller 
problems. 


4. Numerical results 


In this section we conduct a set of numerical experiments in which we vary the amount of demand 
information available to the manager when he makes contingent worker utilization decisions and we vary 
the flexibility of the workers. We investigate how the firm’s costs and the optimal staffing levels are affected 
by information and flexibility. 

We assume that the firm contracts contingent labor from an agency that requires notification of a 
particular day’s labor requirement by a fixed deadline. We assume the total amount of work that arrives 
each day is a Poisson random variable with mean 20. The agency’s notification time splits the total amount 
of work into two batches which we assume are independent Poisson random variables with means that sum 
to 20. The larger the mean of Batch 1 the greater the information available to the manager when notifying 
the agency, in the sense that the variance in the total demand given the Batch | arrival is smaller. Similarly, 
the later the notification time the larger the mean of Batch 1. We consider five information scenarios 
with the following (Batch 1, Batch 2) means: (20, 0), (14, 6), (10, 10), (6, 14), and (0,20). We vary flexibility 
by adjusting the contingent worker guarantee G. For each information scenario we consider eight possible 
values for G, 0, .2, .3, .4, .5, .6, .7, and .8. 

For the purposes of these experiments we model the regular worker absenteeism with a binomial 
distribution. For odd stages t, n, ~ Binomial(N, py) is the number of regular workers present in stage 
t. We are assuming that each worker is present with a probability py, independently of all other 
workers. 

In all the cases we assume the following parameters values: V = 20, Cry = 1, Coy = 1-2, Cot = 2.0, cp = 0, 
Cy = 2.5, cg = 2.5, OT ey = .25, OT wy = .25, Tow = Ror = 1.0, py = .95. 

In Fig. 2 we plot the firm’s total cost under optimal staffing as a function of the contingent worker 
guarantee for the five different information scenarios. We can see that cost increases with decreased flex- 
ibility (as G increases) and costs decrease as more demand information is available earlier. We note that the 
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Fig. 2. Total firm costs for different information scenarios and flexibility levels. 
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Fig. 3. Contingent labor staffing levels for different information scenarios and flexibility levels. 
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Fig. 4. Contingent fraction of workforce for different information scenarios and flexibility levels. 


cost to the firm of staffing with only regular workers is approximately 517 and this places an upper bound 
on costs with contingent workers available. Interestingly we see that flexibility and information each en- 
hance the cost reductions brought about by the other. 

In Figs. 3 and 4 we plot the number of contingent workers contracted and the fraction of the total 
workforce that is contingent, respectively, as a function of the contingent worker guarantee for the five 
different information scenarios. Looking at these two figures together we can see that as information in- 
creases the firm contracts more contingent workers and they make up a greater fraction of the workforce 
for all levels of flexibility. This is interesting because it means that the firm reserves more of the flexible 
capacity as it has more access to demand information. In other words information and capacity are acting 
as complements. This result makes sense because additional information is making each unit of flexible 
capacity more effective. We also see in Fig. 3 that as flexibility increases (i.e., G decreases) the optimal 
contingent workforce size tends to increase. * We note that even when contingent worker decisions are 
made and no demand information is available, the 0—20 scenario, it is optimal to have some contingent 
workers to work down backlog. 

In Fig. 5 we plot the number of regular workers staffed as a function of the contingent worker guarantee 
for the five information scenarios. We see that as flexibility of the contingent worker decreases (i.e., G 
increases) regular full time staffing may decrease. While it is true that the largest full time staff occurs when 
no contingent workers are used we see that increased flexibility of the contingent workers per se does not 
reduce full time employment. This phenomenon is surprising given that it is common to view the increased 
use of contingent workers as a threat to full time employment. Many conclude that anything that would 
make the contingent workers more attractive, like increased flexibility, could only make the negative impact 
on full time employment stronger. In Fig. 5 we see evidence that this is not true in general. The reason 
behind this is that increasing G has two effects. The first effect is that it reduces the number of contingent 
workers you hire, which may increase NV. The second effect is that you want to use the contingent workers 


? We believe that the non-monotonicity displayed in all the graphs is due to integrality effects. 


EJ. Pinker, R.C. Larson / European Journal of Operational Research 144 (2003) 39-55 49 


25 ae | Arrival Pattern: 
4" 7 ond 


Optimal N =o 90-0 
15 aa 14.6 

10 10-10 

5 14-6 

0-20 


0 0.2 0.4 0.6 0.8 1 


Guarantee G 


Fig. 5. Regular worker staffing as function of guarantee for different information scenarios. 
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Fig. 6. Units of backlog carried as a function of guarantee for differing information scenarios. 
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Fig. 7. Number of overtime shifts used as a function of guarantee for differing information scenarios. 
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Fig. 8. Number of contingent worker shifts used as a function of guarantee for differing information scenarios. 


you have hired more since they have been guaranteed work and this may reduce N. When there is a lot of 
demand information, e.g., the 20-0 scenario, the second effect is stronger than the first. 

In Figs. 6-8 we plot the expected total backlog experienced by the system, the expected total number of 
overtime shifts used, and the expected total number of contingent worker shifts used, respectively. We see 
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that backlog tends to increase with less information and less flexibility. This means that both flexibility and 
information make it easier to handle demand variability and the result is lower costs as is seen in Fig. 2. 
Surprisingly, we see in Figs. 7 and 8 that overtime and contingent worker usage are in one sense substitutes 
and in the other complements. When information increases, i.e., moving from scenario 0—20 to scenario 20— 
0, the number of overtime shifts used declines while the number of contingent workers used increases. 
However, as flexibility decreases and the guarantee, G, increases, both overtime and contingent worker 
usage seem to move together. This phenomenon is a result of the way staffing responds to changes in 
flexibility as was seen in Figs. 3-5. 

We have presented results for only one set of cost and productivity parameters. We have conducted 
experiments in which we vary these parameters and generated predictable results. Increasing the cost or 
decreasing the productivity of a worker type reduces the staffing of that type and increases costs. Lower 
backlog penalties lead to higher backlog levels. 

The above example allows us to consider the computing and data requirements for an effective imple- 
mentation of the model presented in this paper. Overall we do not find the model to have data needs beyond 
those managers typically deal with. In theory, the evaluation of different staffing arrangements with the 
model can be computationally intensive since it relies on a dynamic programming model. However, by 
taking advantage of the results discussed in Section 3 it is possible to solve realistic versions of Problem (P) 
with a desktop PC in a practical amount of time. For example determining the optimal staffing level for a 
single information scenario and single value of the flexibility parameter took 70 seconds on a Pentium III 
with a 650 MHz CPU. Experimenting with the performance for larger-scale problems we found that every 
doubling of the system size in terms of expected demand and staff sizes led to an approximately 16-fold 
increase in computation time. The relationship between running time and the parameters V and G is ap- 
proximately linear. Selecting optimal staffing levels would be a planning problem that is typically solved 
infrequently and therefore running times of several hours may be acceptable. On a tactical level computing 
the optimal contingent worker and overtime usage each period for a given staffing level takes only of the 
order of minutes even for problems several time larger than the examples we have presented. 

The parameters of the contingent worker contract: G,OT,, and the regular worker parameter OT, are 
set by the decision maker. The labor cost parameters Cot, Cow, Crw, and cr are also set by the decision maker 
or directly observed. The productivity rates 2.,, 2, although not always accurately determined are 
quantities most managers are comfortable with estimating. The most difficult parameter to determine is the 
penalty for backlog cy. This parameter can have two components to it, a concrete operational cost and a 
less-defined cost as perceived by the firm for providing poor service. In practice we believe that a manager 
should ask herself how much the firm should pay to prevent a unit of backlog with certainty, i.e., eliminate 
the service delay of a group of customers. This assessment should be a reasonable approximation of the 
penalty for backlog. Since it is straightforward to use the model to calculate expected backlog cy can also be 
viewed as a parameter that is adjusted until the system operates with an acceptable level of backlog. The 
final period backlog penalty, cg, can be calculated in a similar manner taking into account any special end 
of planning horizon concerns about backlog. The demand process can be estimated using historical data, as 
can estimates of the absenteeism probability (1 — py). Any changes in the notification times for contingent 
workers would typically require a new analysis of the demand process unless the demand process were a 
stationary Poisson process. 


5. Conclusions 
The management of human resources is one of the most complicated managerial tasks and it is im- 


possible to include all issues relevant to staffing decisions in a single model. We have not modeled turnover, 
training costs or labor supply all of which affect the long term success of any staffing strategy. We have 
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however isolated a primary motivation for utilizing contingent labor, namely demand uncertainty, and 
developed a model that can quantify many of the operational costs of utilizing contingent labor under a 
wide range of circumstances. By using our model a manager can determine the expected short-term cost 
benefits of working time flexibility derived by better matching labor supply to demand. Broader organi- 
zational issues related to contingent labor can be considered in light of these cost benefits. 

We have seen that including the timing of demand information in the model yields a better under- 
standing of the value and impact of labor flexibility. In the numerical experiments we found that infor- 
mation and flexibility act as complements. This means that increasing labor flexibility will not yield better 
performance unless the appropriate information systems are in place to gather sufficient demand infor- 
mation. 

One way to interpret the information scenarios is as different amounts of warning time given to con- 
tingent workers about their work hours for the next period. In the minimal information case workers are 
given the most warning time. Full information gives workers the least advance warning. The numerical 
examples show how the model can be used to quantify the value of work load information or equivalently 
the value of increased flexibility of workers with regard to notification of work hours. As a result, the model 
can potentially serve as a useful tool to resolve, in an equitable manner, the inherent tensions in employers 
and employees needs vis-a-vis working time flexibility. 

Finally we have also shown that the model can be used to solve non-trivial problems with reasonable 
data and computation requirements. This tractability enables the manager to assess the performance of a 
variety of potential contingent labor contracts efficiently and to make quick decisions about their use from 
period to period. 


Appendix A. Proofs 


To simplify the exposition we list a few elementary results on convexity that we make use of in the 
following proofs: 


Property A.1 (Theorem 5.3 in Rockafellar, 1970). If f(x,u) is jointly convex in (x,u) and (x,u) are in a 
convex set S, then min, f(x,u) is convex in u. 


Property A.2 (Theorem 5.1 in Rockafellar, 1970). If f(x) is a non-decreasing convex function and g(x) is a 
convex function, then f(g(x)) is a convex function of x. 


Property A.3 (Theorem 5.5 in Rockafellar, 1970). The positive part ([-]") operation preserves convexity. 


Property A.4. If f(x,y) is convex in x and y, then the expectation of f(x,y) over y is convex in x. This fact 
follows from the fact that multiplication by a positive scalar and summation preserve convexity. 


Property A.5 (Appendix A.5 in Bertsekas, 1976). If f(x) is convex, then the expectation of f (x + b) over b is 


convex in x. 
Lemma A.1. foy(x27, Koy) is a convex function of (x2v, Koy). 


Proof. The final stage, 2V is an even stage in which we make an overtime decision. When cg < Cot /Zot, WE 
never utilize overtime 


=> fry, Koy) = CgxX2y a convex function. 
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When cot /Tot < Cg, We use as much overtime as we can to prevent any backlog 


‘ Xv 
> favor, Koy) = CoMin loris), =| + CR [voy = OT(S)zoi] 
ot 
For x.y < 1OT(S) we have that foy (xo, Kay) is increasing at rate Co./Hor. For xoy > XoOT(S) it is in- 
creasing at rate cg. Since cg > Cot/Mor by assumption, foy is convex. Note also that foy (x27, Koy) is inde- 


pendent of ky. 


Lemma A.2. foy_1(x27_1, Kov_1) is a convex function of (x2v-1,K2v-1). 


Proof. In stage t = 2V — 1, an odd stage, we make contingent worker use decisions 
=> fil%, ti) = MAN { Cew [Ur a Ki)" + Eolfay (fx + 0 - iia K+) |}. 


We can rewrite f,(x,,«,) aS min,,<jo,v) H (x, Kk, u;). The convexity of H() implies the convexity of f,(), by 
Property A.1. To complete the proof of the lemma we need to prove that H() is convex in x,, k;,u,. We can 
write 


A(x, Kt, u;) ca Cow [Ut = Ki] or El fry (A(%, 6, U;), B(Ky, u;))]. 


By Property A.3 both A() and g() are convex in all their arguments. We have seen that fs, does not depend 
on K so we can ignore this dimension. fy) 1s a non-decreasing function in x since more work cannot reduce 
costs and by Lemma A.1 it is convex. Therefore, by Property A.2 foy (A(x, 0, u,)) is convex in (x;, 0, u;). 
Taking expected values with respect to 0 also preserves the convexity of f4y() by Property A.4. Since 
Cow[Uy — K;] is convex, H() is convex in (x;, K;, u;). 


Lemma A.3. foy_2(x2v_2, Koy_2) is a convex function of (x2v~2, Kav-2). 


Proof. When ¢ = 2V — 2, 


Siti, e:) = min {e@, + cy(x; — @:Mot) + Ey, [fist (Xt — @:Tor + O41, Ks)] }- 


MrTot <Xt 
Dropping the stage subscript we can write this cost-to-go function as: f(x,«) = L(x, «) + (Cot/Zot)x, where 


L(x,K) = min {T,(b, «)}, 


[x-OT(S)mot]+ <b <x 


Tot 


T,(b,k) = (« - 5 t+ Ey [frri(d + Pig k)| with b= x — am. 


By Lemma A.2 we know that fi11 (x41, «41) is convex for t= 2V — 2 


=> Ey. [fiai(b+ b41,%)] is convex in (b,x) by Property 5. 


t+ 


Note: k,4; = k, for even stages ¢. Since I’,(b, x) is a sum of convex functions it is convex. 

We can now show that L(x, x) is convex as follows. 

First, define x) = Ayx) + Aox4 for A,, 4. > 0 and 4, + 4, = 1 and similarly Ko = AyKy + Adkd. 

Second, define the following sets: for i € (0, 1,2]: S; = [[x, — OT(S) zo], 2]. 

By the definition of L(): L(x, ko) <I°,(b, «o) for all 6 € So or L(x, Ko) <P (Arb1 + Arbo, Ait) + Ark2) for 
all pairs (b;,b) such that b; € S;. 
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By the convexity of I’,(): 
L(x; Ko) < AT (bi, K}) + AoT (bo, kK) for all pairs (b1, b>) such that Db; € S; 


L(x! <A, min Amin’ ; 
=> L(x, ko) oe (b1,K1) + om (bo, Kp) 


Therefore, L(x, «) is convex, and f;(x;,«,) is a sum of convex functions. 
Lemma A.4. foy_3(x27_3, Kay_3) is a convex function of (x2y_3, K2y_3). 
Proof. In the following we assume that 0 <u,<M. When t = 2V — 3, 


Sil, K) = Min { Cow [dr = Ki” a Fo. fri (De — UjTew + Gail's [Ky = u,)*)]}. 


We can write /;(x,,«,) as follows (dropping the subscripts): f(x, «) = min, {A,(u, «,x)}, where 
A,(u, €,X) = Cow(u — K) + Eolfisi([x — utew + 6]*,0)] ifx<u 
and 
A,(u,«,x) = Eolfisi([x — utew +O)", K—u)| ifk eu. 


We define R1 to be the set of points in the (u, «,x) space such that x < u, and R2 to be the set of points in the 
same space such that x > u. It is easy to show that within each set, R1 or R2, A,(u, K,x) is convex. We prove 
that convexity holds across sets. 

Assume we have two points p; € Rl and p) € R2. For non-negative 4, and /2 such that 2, + 2, = 1 we 
define p; = 2,p, + A2p. such that p; € R1. We first prove that A,(p3) < 2,4,(p1) + 42A:(pr). 

If we define QO; = A,(p3) — 2, A;(p1) — 424,(p2), we have that 


O3 = — Cow(Aaka — Anta) + Bol fess ((Arxr + Aax2 — (Aru + Agu2) tow + 6]*, 0) 
— Ai fii (Pe — Wi Tow 4 0)*,0) Aofisi([X2 — U2Tew + 0)" 


Q3 can be rewritten as 


’ 
> k2 — u2)]. 


O3 = — Cow(Aaka — Agur) + Eol fess ((Arer + Aox2 — (Aru + A2u2) teow + 8)", 0) 
— Ar fess ([1 — U1 Mew + 6)", 0) — Aafert ([x2 — Uottew + OJ", 0) + A2Ay]; 


where 
Ax = fisi([X2 — U2Thew + 0)* ,0) — fis (2 — U2Tlew + lie K2 — U2). 


Since the most that a unit increase in x can reduce cost is coy we can see that AAx < Cew(Aak2 — Aru). 

Because fxy_2 is convex and an increasing function of x (more work in system increases costs) we can 
apply Property A.2 to show that fsy_2([x — umtew + 0]",0) is convex in (x,u). Since 224, < Cow(A2kK2 — Aout) 
we then have that Q; <0. 

We now define ps = 21p; + Ap. such that py € R2, and prove that A;,(ps) < 4, A,(p1) + A24:(p2). If we 
define O4 = A, (ps) = AA; (pi) = AA, (pr) we have that 


On = —CowAi (uy = K1) + Eel fir (Axi + Anx2 = (Aim + AgU2) Tew + Ol, (Aiki + Ark = Ayu, om, AnU2)) 
— Arf ([x1 — Ur Mew +O)", 0) — Aafiss (Leo — U2Tew + 0)", 2 — u2)). 


For similar reasons as Q3, we have that Q, <0. This proves that A,(u,«,x) is convex. This implies that 
A,(u,«,x) is jointly convex in (u,«,x) and therefore we can write: 
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Si(xe, Kr) = min,, A(x, k:,u,), where H(u,,x;,;) is jointly convex in (w,x;,«,). The convexity of H() 
implies the convexity of f,(), by Property A.1. 


Proof of Proposition 1. The proofs of Lemmas A.3 and A.4 only relied on the convexity of the cost-to-go 
functions in the following periods. These were established in Lemmas A.1 and A.2. Therefore we can 
apply the above arguments inductively to show that: f;(x,,«,) is convex, for all t. Note in each of the 
Lemmas A.1—A.4 we have also shown that /;(x;, «,) is the minimum of a convex function of u, for odd ¢ and 
wo, for even t. 


Proof of Proposition 2. {,(«,) is defined as the backlog 5, that minimizes the function I,(b,, «,) defined as 
follows: 
Cot 


D,(b,, Kr) = (4 = =a a5 Ey, Livi (B+ bis Ket)], 


ot 


where b, = x; — Mot @;. 
Assume that £,(0) > B,(x,). By definition I’,(B,(0), «,) — T',(B(i), K,) S 0. Since f.1 (x, K;) is increasing in 
x; (more work implies higher costs) we have that 


We have assumed that cy < Cot/Zor, Which implies that 


B,(«,) — B,(0) > 0 a contradiction. 


The bound on @* follows from Lemma 1. 


Proof of Proposition 3. The result of Proposition 1 implies that when x, = 0 choosing the optimal u,; in odd 
stages ¢, i.e., minimizing f;(x,,0,u,) for u, € [0,M], is a convex minimization problem, where 


Aili, 0, up) = Cowl + Eo,,1 [fer ([X¢ — Uy Tew + Beal 0)] 
Xt 


x 
a (« 7 | + Eo,., fea (te — Uy Tew + Oril*, 0) eg 


Tow Tow 


‘ F — Cow x 
> min, (x,, 0, ur) = min{ X + Eo, f+ ((X a a.".0)} Cg — 


ut cw cw 


where X, = x, — u;Mew and X, € |x, — May, x;]. The statement of the proposition follows. 
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